UCB-PTH-06/15; LBNL-61344 



Landscape Predictions for the Higgs Boson and Top Quark Masses 



Brian Feldstein, Lawrence J. Hall and Taizan Watarf] 

Department of Physics and Lawrence Berkeley National Laboratory, 
University of California, Berkeley, CA 94720, USA 



Abstract 

If the Standard Model is valid up to scales near the Planck mass, and if the cosmologi- 
cal constant and Higgs mass parameters scan on a landscape of vacua, it is well-known that 
the observed orders of magnitude of these quantities can be understood from environmen- 
tal selection for large scale structure and atoms. If in addition the Higgs quartic coupling 
scans, with a probability distribution peaked at low values, environmental selection for a 
phase having a scale of electroweak symmetry breaking much less than the Planck scale 
leads to a most probable Higgs mass of 106 ± 6 GeV for nit = 171 ± 2 GeV. While fluc- 
tuations below this are negligible, the upward fluctuation is 25 jp GeV, where p measures 
the strength of the peaking of the a priori distribution of the quartic coupling. There is 
an additional ±6 GeV uncertainty from calculable higher loop effects, and also sensitivity 
to the experimental value of a s . If the top Yukawa coupling also scans, the most probable 
top quark mass is predicted to lie in the range (174 — 178) GeV, providing the standard 
model is valid to at least 10 17 GeV, with an additional uncertainty of ±3 GeV from higher 
loops. The downward fluctuation is 35 GeV/y/p, suggesting that p is sufficiently large to 
give a very precise Higgs mass prediction. While a high reheat temperature after inflation 
could raise the most probable value of the Higgs mass to 118 GeV, maintaining the suc- 
cessful top prediction suggests that reheating is limited to about 10 8 GeV, and that the 
most probable value of the Higgs mass remains at 106 GeV. If all Yukawa couplings scan, 
then the e, u, d and t masses are understood to be outliers having extreme values induced 
by the pressures of strong environmental selection, while the s, \x, c, b, t Yukawa couplings 
span only two orders of magnitude, reflecting an a priori distribution peaked around 10~ 3 . 
An interesting extension to neutrino masses and leptogenesis follows if right-handed neu- 
trino masses scan, with a preference for larger values, and if Tr and T max scan with mild 
distributions. The broad order of magnitude of the light neutrino masses and the baryon 
asymmetry are correctly predicted, while the right-handed neutrino masses, the reheat 
temperature and the maximum temperature are all predicted to be of order 10 8 -10 9 GeV. 
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1 Introduction 



The Standard Model (SM) of particle physics is both extremely successful and highly predictive. 
It has passed successive hurdles, from the discovery of weak neutral currents to the precision 
electroweak data from 10 7 Z decays. It even explains why protons are so stable and why 
neutrinos are so light. Despite these successes, it is generally not viewed as a fundamental 
theory, but as an effective field theory valid on scales less than about a TeV. This is because 
the Higgs boson mass parameter receives radiative corrections that are quadratically divergent, 
and therefore proportional to A SM , where Asm is the maximum mass scale that the theory 
describes. For large values of Asm, tree-level and radiative contributions to the Higgs mass 
parameter must cancel to a fractional precision of 



where m# is the physical Higgs boson mass. For the SM to be valid up to 5 TeV, a cancellation 
by two orders of magnitude is already required, and to reach the Planck scale requires an ad- 
justment finely tuned to 32 orders of magnitude. Theories that solve this naturalness problem, 
including technicolor, supersymmetry, composite Higgs bosons and extra spatial dimensions, 
have almost defined physics beyond the SM, and are the main focus of potential discoveries at 
the Large Hadron Collider. 

Rather than just being a description of interactions beneath the TeV scale, if the SM is 
valid up to very large energies, for example to the Planck scale Mpi ~ 2.4 x 10 18 GeV, then it 
predicts the mass of the Higgs boson to be in the range 



The lower limit arises from stability of the SM vacuum [TJ [21 [3j H] . A light Higgs boson results 
in the quartic Higgs self interaction becoming negative at large energies; if it is too negative, 
then the universe undergoes quantum tunneling to a phase quite unlike the observed phase. 
On the other hand the upper bound results from perturbativity: if the Higgs mass is too heavy, 
the quartic self coupling becomes non-perturbative at scales well below the Planck scale. It is 
remarkable that the prediction (j2j) receives so little attention, since it closely coincides with the 
current experimental 95% C.L. range 




(1) 



106 GeV <m H < 180 GeV. 



(2) 



114 GeV <m H < 175 GeV 



(3) 
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from direct searches and precision electroweak data, respectively [5]. For a natural SM with 
Asm ~ TeV, the Higgs quartic coupling could lie in a range spanning two orders of magnitude, 
0.03 — 3, whereas for a large value of Asm the range is narrowed to a factor of 3, and consequently 
rriH is constrained to a factor of about \f?>. 

If the LHC discovers physics beyond the SM, or if it discovers a heavy Higgs boson, then 
the apparent agreement between ([2]) and ([3]) will be seen to be a coincidence. However, if the 
LHC discovers a SM Higgs boson in the range of (J2J), and no physics beyond the SM, then 
the case for a very large Asm will be considerably strengthened. In this case the LHC would 
have verified and extended the Little Hierarchy Problem already visible at LEP [B], and there 
would be two obvious interpretations: new physics could appear at several TeV and the SM 
could be accidentally unnatural with apparent fine-tunings at the percent level; alternatively 
the conventional understanding of naturalness could be completely wrong and Asm could be 
extraordinarily large. The LHC would eventually yield a very precise measurement of the Higgs 
boson mass; could this discriminate between these two interpretations? 

An alternative to naturalness is the idea that the universe contains many patches with dif- 
fering underlying physics, and only those patches with a certain complexity will be the subject 
of observation. Thus certain observations can be explained from the selection of complexity 
rather than from symmetry principles [7J. If the weak scale is selected anthropically, i.e. by 
environmental requirements for complexity, the concept of naturalness is not needed to under- 
stand the hierarchy between weak and Planck scales [8]. In this picture, the fundamental theory 
of nature may contain a huge number of vacua, and the Higgs mass parameter may depend on 
the vacuum. If the universe contains many patches, each with its own vacuum, then essentially 
all possible values of the weak scale are realized somewhere in the universe. It is only in those 
patches where the weak scale is in the range of 100s of GeV that atomic physics provides the 
building blocks for carbon based life. The cut-off scale of the SM, Asm, may be extraordinarily 
large, and the Higgs expectation value of most vacua may be of order of Asm, yet such vacua 
are irrelevant to us, since we necessarily find ourselves in one of the very rare patches having 
hospitable chemistry. 

Despite resistance from some physicists, these ideas — a landscape of vacua and environmen- 
tal selection — have gained attention over the last decade: selection of patches of the universe 
containing large scale structures, such as galaxies, can essentially solve the cosmological con- 
stant problem and predicts that observers are likely to inhabit patches containing dark energy 
[9]. Furthermore, the number of known vacua of string theory continues to increase, and the 
string landscape appears to be able to scan parameters sufficiently densely to allow environ- 
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mental selection [TU|. Such ideas should perhaps be resisted, since they are so difficult to test 
experimentally. How are we to test such theories when the other patches of the universe lie 
outside our horizon? Natural theories are tested by assuming a theory with a parameter space 
that is sufficiently restricted that predictions for observables can be made. In the landscape, 
predictions may be possible by combining environmental selection with assumptions about the 
underlying vacuum probability distributions. In this paper we argue that the Higgs boson mass 
is well suited to this sort of prediction, since it depends on the quartic scalar coupling, and 
if this coupling is too low the SM vacuum decays. If the landscape favors low values of this 
coupling, environmental selection will favor Higgs masses close to the stability bound. 

Expanding the Higgs potential V(<j>) in powers of (ft/ Asm, the relevant part of the Higgs 
potential is 



In this paper we assume that the entire SM Higgs potential varies, or scans, from one patch of 
the universe to another. The scanning of A 4 leads to the scanning of the cosmological constant 
and allows an environmental selection for galaxies, largely explaining the observed value of the 
dark energy [9] . Similarly, the scanning of the Higgs mass parameter m 2 and quartic coupling A 
allows for the environmental selection of required atomic properties, determining the weak scale 
a/— m 2 / A [8] . In those regions of the universe selected environmentally by both the cosmological 
constant and the weak scale, there will be a variation in the coupling A and therefore in the 
Higgs boson mass from one patch to another. What are the physical consequences of this 
scanning of the Higgs quartic coupling A? Environmental selection will choose the required 
electroweak phase: we must live in a region with A above a minimum value, A c . This is how 
the usual stability limit on the SM Higgs boson mass arises in the landscape, leading to the 
lower bound in ([2]). For A > A c there is no obvious environmental selection of one value of 
iriH over another, hence the prediction for the Higgs boson mass is governed by the probability 
distribution^) P(X) of the coupling A in the vacua with acceptable large scale structure and 
atomic physics. To proceed, apparently one needs a calculation of -P(A) from the landscape. 
In fact, only a single assumption on its form is required: we assume that P(X) is sufficiently 
peaked at low values that A is expected to be near A c , as illustrated in Figure [TJ corresponding 
to a Higgs mass near the stability bound. The precision of this prediction clearly depends on 
how steep -P(A) is near A c , as we investigate in some detail. 

It is important to stress that our prediction does depend on an assumption for P(X). For 

2 To be more precise, by the probability distribution P(m 2 , A)| / _ m2 /A _ 174 GcV - Hereafter, we mean this by 
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P(A). 
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Figure 1: The a priori probability distribution is peaked at low values of A. 



example, an alternative assumption that P(X) is peaked at large positive values of A predicts m# 
near the perturbativity limit of 180 GeV. However, this prediction could depend on whether the 
theory has a large energy interval in which it is strongly coupled before a perturbative quartic 
emerges at low energy. Furthermore, an observation of a Higgs boson mass in this region could 
be interpreted in terms of a strongly coupled fundamental theory without any recourse to the 
landscape. In the rest of this paper, we concentrate on the possibility that the probability 
distribution P(X) is peaked towards low values, corresponding to a Higgs boson mass near the 
stability limit of m# ~ 106 GeV. A measurement of such a Higgs mass, together with the 
absence of any new physics, would provide evidence for environmental selection. 

The top quark Yukawa coupling, h, plays an essential role in electroweak symmetry breaking 
in the SM with large Asm, through its effect on the renormalization group (RG) scaling of the 
Higgs quartic coupling, A. Suppose that h also scans in the landscape, so that there is a 
combined probability distribution P(X, h). How is the discussion of the Higgs mass prediction 
changed? As we will show, the metastability boundary in the (A, h) plane has a special point, 
where A is at a minimum, and a simple assumption for the a priori probability distribution, 
P(X,h), implies that patches of the universe in the desired electroweak phase are most likely 
to be in the neighborhood of this point. Remarkably, at this special point 

m Hc = (121 ± 6) GeV m tc = (176 ± 2) GeV, (5) 

where the uncertainty corresponds to Asm — 10 18±1 GeV, leading simultaneously to a prediction 
for the Higgs mass and to a broadly successful post-diction for the top mass. Although there 



are further uncertainties of about ±6 GeV and ±3 GeV on m# c and m tc respectively, from 
higher loop effects, this result is nevertheless very striking^ 

It is well-known that the stability limit on m# in the SM depends on Tr, the reheat tem- 
perature of the universe after inflation. In the landscape, Tr may scan from one patch of 
the universe to another, or it might be fixed; either way, we study the landscape Higgs mass 
prediction in two cases. In the first case, in sections 2 and 3, we assume that the probability 
distribution is dominated by patches having Tr ^ 10 8 GeV, so that thermal fluctuations can 
be ignored, and the stability limit arises from quantum fluctuations from the false vacuum. In 
section 2 we keep the top Yukawa coupling fixed, while in section 3 we allow it to scan. In the 
second case, in section 4, we assume that the landscape probability distribution is dominated 
by patches having Tr ^ 10 8 GeV, and study the Higgs mass prediction arising from thermal 
nucleation of bubbles of the true vacuum. In section 5 we embed these ideas in the "scanning 
SM" where all SM parameters scan, discussing charged fermion mass hierarchies, the scale 
of neutrino masses, leptogenesis and consequences for inflation. Conclusions are presented in 
section 6. 

2 Higgs Mass Prediction from Quantum Tunneling 
2.1 Quantum Tunneling 

In this section we predict the Higgs boson mass using the environmental constraint of sufficient 
stability of the electroweak vacuum in a landscape scenario. The SM is assumed to be valid 
up to a high energy scale, Asm, such as the Planck scale. To be precise, we define Asm to 
be the scale at which the standard model RG equations and bounce action calculations are 
still valid to within 1%. Thus Asm is slightly lower than the scale at which new physics 
actually arises. We concentrate on the consequences of scanning (A 4 ,m 2 ,A), assuming that 
all other parameters of the SM do not effectively scan or are fixed tightly to observed values 
by other environmental constraints. In the patches of the universe dominating the probability 
distribution, the maximum temperature after inflation is assumed to be sufficiently low (Tr & 
10 8 GeV) that thermal nucleation of the phase transition is sub-dominant. 

At large field values, the effective potential of the Higgs field H = \/2Re(0°) is well approx- 

3 We should mention here that predictions for the Higgs and top masses in the SM, assumed valid to very 
high scales, have been made based on very different principles. One prediction resulted from assuming that 
there are two degenerate electroweak vacua, with one occuring near the Planck scale [llj , while another resulted 
from the assumption that two phases should coexist, giving borderline vacuum metastability |12| . 
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imated by 

V(H) = ^Ih* for H^>v = 246 GeV, (6) 

where A (if) is the value of the running Higgs quartic coupling evaluated at the scale \i = H . 
The RG equation for A has the 1-loop form 

16vr 2 ^ = 24A 2 + I2\h* - 6h* - 9Xgl - 3Xg Y + \g Y + \g\g Y + \gt (7) 

where h is the top Yukawa coupling, and gi and gy are the SU(2)l and U{l)y gauge couplings, 
respectively. Several trajectories for A(/i) are illustrated in Figure [2j The effect of the large top 
Yukawa coupling is to reduce the value of A at high energies. In fact, A will become negative 
below Asm if its value at \i ~ v is too small. In such cases, the minimum with (H) = v is 
not the true minimum of the potential, and we must consider the possibility of the universe 
tunneling quantum mechanically from our hospitable (H) = v vacuum to the inhospitable one 
with very large (#)@ Throughout we assume that conditions of the early universe lead to a 
sufficient number of patches with the desired metastable phase. 

The quantum tunneling rate per unit volume to the true vacuum is dominated by a bounce 
solution H{r). A pure up-side-down quartic potential, ([6]) with negative A, has an S0(4) 
symmetric bounce solution [13]. Due to the conformal nature of the potential (ignoring the 
running of A and the quadratic term of order the weak scale), there is a family of bounce 
solutions, H(r) = cH (cr) for c > 0, with different size oc 1/c and field value at the center 
H(r = 0) oc c. A bounce solution with H(r = 0) = M contributes to the decay rate an amount 

r(M)^M 4 e"^W, (8) 

so that the total decay rate is approximated by [H [21 [31 H] 

(9) 



TfAol ~ max T(M) = max 

M<A 5M M<A SM 



Stv 2 

j\^f 4 e ~3|A(A/)| 



Each patch of the universe has its own A = A(AsAf), RG trajectory A(/i), and corresponding 
decay rate T[Ao]. By time t, an arbitrary point remains in the desired false vacuum with 
(H) = v provided no bubble nucleated in its past light cone. Since the vacuum tunneling rate 
r[A ] does not depend on time, bubble nucleation is most likely to occur at the epoch with time 
of order t, rather than at a much earlier epoch. Assuming that the universe has been matter 



4 In this phase the ratio of quark masses to the QCD scale and to the Planck mass are very far from allowing 
conventional nuclei and cosmology. 
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dominated most of the time until t, the fraction of the volume of a patch in the false vacuum 
at time t is given by [Hj 

/(A ;t) = e- r ^ i4 , (10) 

where a coefficient of order one in the exponent is ignoredM^ 

For practical applications, we are interested in t of order 10 10 years; t is either still in the 
matter-dominated era or in the dark-energy dominated era so that the exponent is roughly of 
order r[A ]T 4 , where r = 10 10 years. If r(M)r 4 ^ 1 for any scale in the range v <C M & Asm, 
the survival fraction f(Xo] t) is significantly less than 1. Since T(M) primarily depends on A(/i) 
renormalized at \x = M, one can introduce a critical value of the quartic coupling A C (M) by 
r(M)r 4 ~ 1 

A C (M) K { - I „-) « -0.047 ( 1 + -Lin (11) 

where ln(Mp^r) « 138 is used to obtain the last expression. This critical line of stability, A C (M), 
is shown by a dashed line in FigureEJ Whenever a trajectory A(/x) dips below A c (yu) at any scale 
v <C ^ A 5M , bounces with i7(r = 0) of order that scale would have destabilized the desired 
false vacuum. The usual metastability bound on the Higgs boson mass in [21 [31 H] corresponds to 
the trajectory that touches the critical line of stability at one scale \x = Md om ~ 10 16 GeV, while 
maintaining X(ji) > A c (/i) for all \i up to Asm- Thus, the bounce of most danger is the one with 
H(r = 0) = Mdom, and the metastability bound is independent of Asm provided Asm > Mdom- 
Note that there is an assumption here: for all trajectories of interest, the contribution to T[Ao] 
from quantum fluctuations to field values larger than Asm are sub-dominant. We must assume 
this because we cannot calculate these contributions without knowing some of the details of the 
more complete theory above Asm- The trajectories of interest include the special one, shown as 
a solid curve in Figure [2j and all those that lie above it. The assumption is reasonable because, 
as the SM begins to break down at energies just above Asm, these trajectories all lie above the 
critical curve A c (/x), shown dashed in Figure El 

A landscape will provide some a priori probability P(Xq), and we follow the principle that 
the probability distribution of observable parameters is further weighted by the fraction of 

5 The exponent depends non-perturbatively on the Higgs quartic coupling A as in ([8]), and, if |A| is small, a 
small variation in the value of A can change T by orders of magnitude. Thus, any prcfactor of order unity to 
Tt A is unimportant compared with the sensitivity on A. 

6 As the universe evolves into an era dominated by the cosmological constant the fraction of the volume in 
the false vacuum is given asymptotically by /(Ao;i) = e~ r ^°' t ' H ° , where Ho is the Hubble parameter and a 
coefficient of order one in the exponent is again ignored. 
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Figure 2: Several RG trajectories for the Higgs quartic coupling, A, are shown by dotted curves. 
The dashed line gives the critical coupling A c at each energy, defined by ( TTTT) . The special tra- 
jectory that corresponds to the usual SM metastability bound is shown by a solid curve. Any 
patch of the universe with a trajectory that passes significantly into the region below the dashed 
line is likely to decay out of the false vacuum with (H) = v before a cosmological time scale 
of order 10 10 years. In this figure, Asm is chosen to be 10 18 GeV, and (gy, gi, g s , h)\^ = \ SM = 
(0.466,0.512,0.500,0.391) yielding, for the special trajectory with A (A<? M ) = -0.0462, weak- 
scale SM parameters (a(m^) -1 , sin 2 9w(mz), o^ s (mz), m t ) = (127.9, 0.2313, 0.1176, 171 GeV) 
|T5] as well as m# = 106 GeV. The analysis is performed with 2-loop RG equations [16] 
and 1-loop (2-loop in g s ) threshold corrections [TTJ HH], [19] at the matching scale fi = m t [20] . 

observers who see them. For patches of the universe having X(Asm) = Ao, only a volume 
fraction f(Xo]t) remains in our desired false vacuum with (H) = v, introducing an additional 
A dependence in the total probability distribution 



The a priori distribution -P(Ao) describes all Ao dependence that originates from cosmology 
up to the end of (the last) inflation. Since the Higgs boson mass itself is not environmentally 
important, / will be the only Ao dependence (and Higgs boson mass dependence) that originates 
after the end of inflation. 

One could also study an all-time probability distribution, instead of the distribution for 
contemporary observers. Such a distribution is obtained by integrating (fT2|) over t, weighted 
by the time evolution of the number of observers, p{t), which we expect to be independent of 
A . Thus, / still remains the only source of the late-time dependence of the distribution on the 
Higgs boson mass. 



V(\ ;t)d\o = P(\o)e- T[Xo]ti 



d\ . 



(12) 



S 



2.2 The Prediction 



Using the probability distribution of (|12p we would like to answer three key questions. 

1. What is the most probable value of Ao, and therefore the Higgs mass, as seen by observers, 
such as us, living at the present age of the universe t = to ~ 1.4 X 10 10 years? 

2. What range of Higgs masses correspond to a probability within a factor of, say, 1/e of the 
peak probability? 

3. What fraction of observers in the "multiverse" live at the time t or later? 

The answer to the first question gives the central value of the Higgs mass prediction, while 
the answer to the second gives a measure of the uncertainty of the prediction. We will be 
particularly interested in how these two answers depend on the form of the a priori probability 
distribution P(A ). We suspect that a precise prediction will follow only if there is some peaking 
in this distribution. There may be a trade-off: if the distribution is too weak, then we cannot 
predict the Higgs mass with significant accuracy; if it is too strong, then we may find ourselves 
living in an unstable patch of the universe that was extremely lucky to survive until now. With 
the third question, we thus investigate how much peaking we can tolerate before observers 
living as late as to become a rarity. 

We assume that Asm ~ M dom . Since tunneling is most likely to occur at the scale M dom , 
as can be seen from figure^ with this assumption the analysis becomes independent of Asm- 
Consider patches of the universe with various values of X(M dom ) as illustrated by the trajectories 
in Figure [2j An a priori probability distribution peaked at low values of Ao, as illustrated in 
Figure [U means that there are more patches with X(M dom ) below the critical value X c (M dom ) 
than above. However, these patches are more likely to decay. The most probable observed 
value of A today is therefore determined by a competition between -P(Ao) and e~ rt ° in (1121) . and 
is given hy^\ 

X(M dom ) = X c (M dom ) ( 1 - In — ) (13) 

\ 54U p J 

where the parameter p describes the strength of the peaking of the a priori probability distri- 
bution at Ao = \ c (Asm) = A Cj o and is defined by 



<91nP(A ) 

p =- 



d In A r 



(14) 

-V>=-V),c 



7 We ignore the difference between to and r = 10 10 years, and use the latter as a reference, rriH in (TT5 
increases no more than 0.1 GeV by using t$ instead (e.g., [3]). 
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Since A c is negative, we are interested in positive values for p. (Note that we use A , c and A Cj0 
interchangeably.) Clearly, except for extremely large values of p, the double exponential form 
of the decay probability factor is so powerful that we can just take A(M c fo m ) = A c (Mrf om ). Thus 
the most probable Higgs mass measured by observers at time to corresponds to the solid line 
trajectory in Figure [2j which also corresponds to the usual SM metastability bound. We find 
that this trajectory corresponds to a Higgs boson mass§ 

fn H = 106 GeV + 6 GeV ( ^~]J 1G ^ ) _ 2 . 6 GeV ( °» ~ °- 1176 > | ± 6 GeV . (15) 

V 2 GeV J V 0.002 J v ; 

There are several theoretical uncertainties that contribute to the ±6 GeV uncertainty quoted 
in (115I) . in particular in relating the physical Higgs boson mass to the running quartic coupling 
evaluated at very high energies and in relating the physical top quark mass to the running 
top quark Yukawa coupling. The present 6 GeV uncertainty in % coming from the ±2 GeV 
experimental uncertainty in m t will be reduced by a factor of 2-4 at the LHC, depending on 
how well systematic uncertainties can be understood, while a reduction in the uncertainty in a s 
may have to wait for the ILC. Finally, while rf%H is the most probable value for us to observe, 
a crucial question is the width of the probability distribution for the Higgs mass, which will 
depend on the strength p of the peaking of the a priori distribution P(Xq). 

Reference [2TJ has studied extensions of the SM where m 2 scans but not A. It was assumed 
that Ao is negative so that the RG trajectory for A passes through zero at some scale A cross . 
Only those patches of the universe with m 2 negative and |m 2 | < A 2 ross will have a metastable 
hospitable electroweak phase, and hence environmental selection will lead to the weak scale 
being close to A cross , giving an alternative understanding of a low weak scale. This is a very 
interesting idea, but we stress that it is very different from the scanning of A that this paper 
is based on. From Figure [2] it is clear that their idea cannot work in the SM, hence in addition 
to a light Higgs boson, they predict other heavy states at the weak scale. 

2.3 The Width of the Prediction 

In the SM the Higgs mass must lie in the experimentally allowed range of (114-175) GeV. To 
see evidence for a landscape, it is crucial that the probability distribution for the Higgs boson 
mass is sharply peaked compared to this allowed range. How strong an assumption does this 
require for the a priori probability distribution? 



8 These results are in agreement with the metastability bound found in [1] . The results in previous versions 
of our paper gave values of mu about 5 GeV larger (for a given m t ) due to incorrectly neglecting part of the 
top threshold in h t . We thank G. Giudice for a very useful e-mail 
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To study this question we Taylor expand the probability "P(A ) oc Pe~ Tri of (TT21) about the 
most probable value of the quartic coupling A , and find the couplings A± where V falls to 
1/e of its maximum value. We choose A_ < A + , so that A + corresponds to the heaviest Higgs 
boson that we are likely to observe. The maximum value of the exponent Tr 4 is p/540, which, 
for reasons discussed in the next sub-section, we take to be less than unity. It follows that the 
curve V{\q) is highly asymmetrical about Ao, with A_ determined by the very rapid drop off 
of e _r? ~ 4 

A-=A„fl + 4ta-V (16) 



and A + determined^! by the drop off of P 



540 p J 



A + = Ao (l - i) . (17) 

We immediately see that, to obtain a significant Higgs mass prediction, there is no need to 
have a very large p; a much more modestly peaked a priori distribution is sufficient. We also 
see that for modest values of p, the uncertainty in the Higgs mass coming from the landscape, 
i.e. from the unknown a priori probability distribution and from the probabilistic nature of 
quantum tunneling, is dominated by the 1/p contribution to A+, which increases the Higgs 
boson mass above m#. The difference between A_ and Ao corresponds to a change in the Higgs 
mass that is negligible compared to the uncertainties of (fl~5l) . Hence, while m# of (fTBT) is the 
most probable value for the Higgs mass, the expected range is from rf%H up to where 
Aran is determined from the shift in the quartic coupling at Asm' AA ~ |Ao, c |/p ~ —.05/p. 
RG running will not change this shift by much as it is taken to the weak scale, and since 
m H ~ 350^ GeV, we find@ 

25GeV 

Am H « . (18) 

p 

A very precise prediction follows if p > 10, while a reasonably precise prediction requires only 
p > 3. In the more complete theory at the scale of new physics, the description of the relevant 
landscape of vacua may involve a probability distribution with some other set of couplings. 
The simple peaking of P, corresponding to p > 3, may have a simple origin in the distributions 
of the more complete theory, because such a peaking is not destroyed by mild RG scaling. 



9 We caution the reader that on Taylor expanding P about Ao, the expansion is on the verge of breaking down 
on reaching A+. We have checked that (fT7|) is precisely accurate for the cases of power law and exponential 
distributions, but in general it should only be taken at the factor of 2 level. 

10 A different expression is required for p <; 0(500), but for such large values, the predicted range is so narrow 
that the precise width is unimportant. 
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2.4 The Fraction of Observers Living After t 

There is a second important consequence of having the a priori probability distribution strongly 
peaked at low values of Ao: it becomes more probable for observers in the multiverse to live soon 
after non-linearities in the large scale structure appear. Once again, this can be understood 
by studying figure [2j As lower values of Ao become more strongly preferred, so more patches 
of the multiverse have X(M dom ) < X c (M dom ), and these patches typically decay before time t . 
This is clearly a problem if there are so few patches of the universe that survive until to, that 
we are very rare observers. How does this limit the degree of peaking of P(Ao)? 

To get a feel for this we can study how P varies as A is decreased below A c . Instead of 
using Ao as the variable describing each patch we can use t, the time at which the patch with 
coupling A typically decays. Taylor expanding about t to first order, we find 

P(t) =P(*o) (l + ^ m (19) 

The earliest relevant time is the time that suitable structures first went non- linear, tNL, and we 
expect that ln^o/tAM.) is n o more than a few. Hence we see that observers at t are common 
in the multiverse providing the peaking of the probability distribution satisfies 

140 , , 

V < j 77 77 y 20 

\n(t /t NL ) 

which certainly allows for a precise Higgs mass prediction, as can be seen from (TT8]) . 

Even for larger values of p the probability for observers in the multiverse at time to does 
not drop extremely quickly. Approximating the time evolution of the number of observers in 
an undecayed patch, p(t), as a step function at t NL , the fraction of observers in the multiverse 
living at time to or later is 

(21) 

where 

POO 

I(t) = / P(A )e- r[A(,] *' 4 dt'd\ . (22) 



To evaluate this expression an explicit form for P is required. For example, for an exponential 
distribution 

P = Ce p ^ (23) 

we find 



f w i if ' (24) 
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so that the fraction of observers at to is n °t highly suppressed even for p as large as several 
hundred. For a power law distribution, the times t in ( l24l) are replaced with InAsMt, allowing 
considerably larger values of p. To summarize: there is a wide range of p that gives a precise 
Higgs mass prediction while allowing observers at to to be common. 

3 Scanning the Top Quark Yukawa Coupling 

In the previous section we made an important assumption: the only SM parameters scanning 
in the landscape were those in the Higgs potential. Of all the other SM parameters, the one 
in our patch of the universe that has the most importance for electroweak symmetry breaking 
is the top quark Yukawa coupling, h, via the RG equation for the scalar quartic coupling ([7j). 
Hence we now discuss the consequences of allowing h to scan, leaving a discussion of allowing 
other Yukawa couplings to scan until section 5. 

Remarkably, we will show in this section that, allowing h to scan, we can go further than 
merely preserving the Higgs mass prediction of section 2; using a landscape probability distri- 
bution, -P(Ao, ho), peaked in Ao beyond the SM metastability boundary, we will be able to add 
a successful understanding of the top quark mass to our prediction for the Higgs mass. 

In this section we again assume that the temperature was never above about 10 8 GeV, 
at least for the patches of interest that dominate the probability distribution and evolved to 
the desired metastable phase, so that we may ignore effects of vacuum transition by thermal 
excitation. 

3.1 The Metastability Boundary in the X—h Plane 

When only Ao scans, the metastability boundary is given by the critical value of Ao that cor- 
responds to \(Mdom) = A c (Md om ). If Ao and h both scan, what is the critical line for the 
metastability boundary in the \o-ho plane? As noted in the introduction, it has a special 
shape that will make an understanding of the top quark mass possible. 

Consider a scale A < Asm, and fix a value of A(A). Figure Eh shows trajectories of A, 
starting from this initial condition, for various values of h. In the figure, A = M dom was used, 
but the qualitative picture is independent of this choice. When the top Yukawa is too small 
or too large, A runs into the forbidden unstable region. As A(A) becomes smaller, the allowed 
range of h shrinks accordingly. When A (A) actually lies on the metastability boundary itself, as 
in Figure[3b, the only permissible value of h remaining is the one for which A runs parallel to the 
boundary at the chosen scale, trajectory c of Figure [3b. The resulting metastability boundary 
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(a): A (A) = -0.0477 

A 




14 15 16 17 18 

X (b): A(A) = -0.0487 

-0.044 r I 




Figure 3: RG trajectories for A(/z), for different values of h and a fixed value for A(A), shown 
for a few decades in energy below Asm — 10 18 GeV. The scale A is chosen to be 7.4 x 10 15 
GeV, corresponding to the value of M^ om found in section 2 for the special trajectory with 
Ao(A,sm) = —0.0462 that gives m t = 171 GeV and m H = 106 GeV. The trajectories in the upper 
panel labeled A, B, C, D and E have A(A) = -0.0477 with h(A) = 0.470, 0.447, 0.419, 0.391 and 
0.367, respectively. Those in the lower panel, labeled b, c and d, are for A(A) = A C (A) = —0.0487, 
with h(A) = 0.447, 0.421 and 0.391. In both panels the dashed line is the line of metastability. 
The analysis uses SM RG equations at the 2- loop level. The gauge couplings at Asm are n °t 
scanned, but are fixed to the values (gy, g^, g s )(AsM) = (0.466, 0.512, 0.500) for all trajectories, 
the same as in Figure [21 Uncertainties from higher loops and from the experimental uncertainty 
on a s are not included. 
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in the X(A)-h(A) plane is shown in Figure Hh- The trajectories labeled A, B, C, D, E and b, c, d 
in Figure [3] appear as points in Figure HI The special trajectory c in Figure [3b is clearly at a 
special location on the met ast ability boundary in Figure UK If there is a probability pressure in 
the landscape pushing A (A) towards A C (A), as in section 2, then we may simultaneously obtain 
the prediction that h(A) will lie at the critical value corresponding to the trajectory c. 

Now, the low energy value of the top mass corresponding to the critical point depends on 
the actual scale A we choose. If A is increased, the general shape of the metastability boundary 
remains the same, as shown in Figure Hb for A = 10 18 GeV, but the special trajectory has 
changed from c to near B, with a consequent change in the top mass prediction. The prediction 
for h should be made at the scale at which there is peaking in the probability distribution 
pushing towards the critical point. Although with h scanning, this is not a RG invariant 
requirement, it is a reasonable possibility for this to occur close to the scale of new physics, in 
particular, at the SM cutoff Asm- What is remarkable is that if we take A = Asm to be within 
a few orders of magnitude of the Planck scale, then the corresponding prediction for h yields a 
value for the the top mass that agrees closely with experiment. The uncertainty from varying 
A — Asm is mild. 

The astute reader will notice a problem with this analysis; the A — h metastability boundary 
was derived using the SM RG equations and tunneling rate calculations, which are expected 
to break down above the scale Asm- As a result of effects from the new physics, the Xo-ho 
metastability boundary will differ from the form shown in Figure HI and is sketched in Figure [51 
However, the possible deviations are greatly restricted: only the half of the curve to the right 
of the critical point, corresponding to tunneling events to scales beyond Asm, will be affected. 
Moreover, this part of the curve cannot change in an arbitrary way; Ao, c is still a strict lower 
bound on the Higgs quartic at Asm, an d thus the true curve must lie above the dashed line of 
Figure [5j As long as the true curve rises somewhat above this horizontal line, the central value 
of our prediction will be unchanged. If this curve did remain too flat, however, our prediction 
would have a large error bar on the high side; m t would have a lower bound, but not much of an 
upper bound. This situation could arise if tunneling events to scales larger than Asm are very 
suppressed!^] To obtain a well controlled prediction for m t , we will thus make one additional 
assumption in this section: either 

11 The tunneling rate depends exponentially on the action. The reason that the metastability boundary 
changes so gradually in the SM is because of tree-level scale independence. However, the tunneling rate of the 
full theory is likely to differ substantially from that of the SM, so that we expect the actual behavior to be more 
extreme than the two dashed curves of Figure [5] 
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(a) A = M dom = 7.4 x 10 15 GeV 
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(b) A = Asm = 10 18 GeV 
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Figure 4: The metastability boundary in the A(A)-/i(A) plane for (a) A = M^ om , found for 
the metastability trajectory of section 2, and (b) A = Asm = 10 18 GeV. The points labeled 
A, B, C, D, E and b, c, d correspond to the trajectories of Figure [31 The metastability boundary 
is shown dashed in regions where the calculations unreliable; using SM equations, the dominant 
scale for tunneling is found to be above Asm = 10 18 GeV. 



1. The metastability boundary rises sufficiently rapidly above h 0)C so as to provide a reason- 
ably sharp upper bound on the top mass, or 

2. There is a fundamental probability distribution pressure in the landscape pushing ho 
towards smaller values, and hence towards h 0>c . We take q = — 9ln Qi^^ h °^ U an d require 
1 ^ q < p. 



3.2 The Critical Value of the Top Mass 

The critical value of the quartic coupling at any scale /x, A c (/i), is defined by r(/i)r 4 = 1, where 
t = 10 10 years. By taking a derivative with respect to \i and using (jSJ), we obtain 

16tt 2 -^- = 2A\ 2 C . (25) 
am/i 

On the other hand, the critical RG trajectory for A, with h scanning, touches the curve A c (/z) 
tangentially at Asm- In this way, the condition for the critical top Yukawa coupling h c follows 
from 

(26) 



c/ln/i 

giving 



\ -\ ,,-a din a 



M=A< 



(\2\ c hl - 6hj - 9\ c gl - 3X c g 2 Y + ~g* + \g\g\ + -g\\ = 0. (27) 

> ' Asa/ 

Since A Cj o < 0, the first 2 terms in (|27|) are negative while all the rest are positive. Taking 
A Cj0 from equation (ITTj) . and using two loop SM RG running for QiXS-sm) an d ^(^54/)^ we 
find the critical value for the top coupling h C) o = Ii c (Asm)- For example, for Asm = 10 18 GeV, 
h c ,o — 0.417. RG scaling to the weak scale, we find a critical top quark mass 



A 



m tc = (176.2 ±3 + 2.2 log 10 ) GeV. (28) 

where the uncertainty ±3 GeV comes from higher order contributions to RG scaling and to 
threshold corrections at the scale of the top quark mass. This critical point is denoted by an 
open circle in Figure [5j Contours of m t and mn are also shown in this figure by the near 
vertical blue dotted and red solid lines. Clearly a crucial question is how far from the critical 
point a typical patch in the false vacuum is likely to be. 



12 The choice of the Higgs boson mass as a boundary condition of the RG equation has negligible impact on 
9l{^sm) and gy(A-SM)- 
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h>o,c 

Figure 5: The metastability boundary in the A -/io plane for Asm — 10 18 GeV. For h < h 0tC 
the solid curve is reliably computed in the SM. For h > h c , the metastability boundary must 
lie above the dashed horizontal line. Two possible boundaries are shown; the lower (upper) 
dot-dash curve results when tunneling in the full theory is suppressed (rapid). The nearly 
vertical blue dotted lines are contours of the top mass in 2.5 GeV intervals, while the red solid 
lines, also not far from vertical, are contours of the Higgs mass with 10 GeV intervals, shown 
for a s = 0.1176. The open circle is the critical point, giving the most probable observed values 
of m t and m# in the multiverse, corresponding to trajectory B of Figure Hb. The closed circle 
represents a patch of the multiverse on the metastability boundary where the top quark is 
measured to be 171 GeV. In our patch, the top mass is within a couple of GeV of this number, 
and hence lies somewhere between the m t of 170 GeV and 175 GeV contours. The importance 
for the Higgs mass prediction of decreasing the experimental uncertainties on m t is clear. An 
arrow denotes the magnitude of the expected upward fluctuations of the quartic coupling A 
above its critical value for peaking parameter p = 20. Such a fluctuation above Ao, c would lead 
to a point on the metastability boundary where the top mass is about 171 GeV and the Higgs 
mass about 110 GeV. 
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3.3 The Width of the Top Mass Prediction 

We would next like to determine how far from m tc the top mass could go, Am t ±, before the 
probability would have fallen by a factor of e. In order to lower the top Yukawa ho below its 
critical value ho >c , the Higgs quartic must be raised to prevent the RG flow of A from heading into 
the forbidden region. This costs probability through the Higgs quartic probability distribution. 
As in section 2, we lose a factor of e in probability when the Higgs quartic is raised above the 
metastability boundary by an amount \X c \/p. The lower bound on the top mass prediction then 
corresponds to the smallest top Yukawa coupling compatible with this raised value of A, and 
can be read off from the solid curve of Figure For a small deviation from the critical point 
(Ao,c, h 0jC ) = (—0.047,0.417), the metastability boundary rises as (A — A 0|C ) — l.2(h Q — h 0tC ) 2 , 
and the downward shift Ah ~ y/\\ 0tC \/(1.2p) or A/i 0)C //i 0)C ~ 0.47/^/p gives the lower bound. 
Running down to the weak scale, Ah/h shrinks by a factor of 0.4, and we find a corresponding 
downward shift in the top mass of 

35GeV , . 

Am 4 „ « — . 29 

VP 

Such a shift is illustrated and discussed in Figure [5] for p = 20. 

The upper bound on the top mass prediction is less certain; it depends on either the shape 
of the A — h metastability boundary to the right of the critical point, or on the strength, given 
by q, of the the fundamental probability distribution on h. In the former case, the upper bound 
depends on unknown physics beyond Asm and we cannot make a quantitative statement. In 
the latter case the probability will have dropped by a factor of e for a top Yukawa that is above 
the critical value by an amount ~ h c j q, corresponding to 

75 GeV 

Am t+ « . (30) 

Q 

Putting everything together, our final result for the top mass prediction is 

m < = ( 176 2 ± 3 + 2 2 log - ife - -j=m + ps) GeV - (31) 

We find this to be a striking success. Indeed, for the agreement with experiment not to be 
accidental, the a priori distribution for Ao must be highly peaked, certainly with p <; 20, so 
that the uncertainty in the Higgs mass prediction from the landscape, ffl8|) . is less than 1 GeV. 
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4 Higgs Mass Prediction from Thermal Fluctuations 



All the analysis in the preceding sections assumed that the temperature of the thermal plasma 
remained sufficiently small after the last era of inflation. On the other hand, thermal fluc- 
tuations can themselves induce vacuum transitions from the false vacuum with (H) = v to 
the true vacuum with large (H). If the temperature was sufficiently high, then the vacuum 
transitions are dominated by these thermal processes, not by the quantum tunneling discussed 
so far [22l [23l [18], [21] . This section investigates how the results in the earlier sections should be 
modified in this case. We will continue to assume that the Standard Model is a good approxi- 
mation up to a very high scale, and will return to the assumption that the top quark Yukawa 
coupling does not effectively scan in the landscape. 

4.1 Vacuum Transitions at High Temperatures 

As the inflaton decays, the inflation vacuum energy density Vi is converted into a thermal 
plasma. The maximum temperature of this plasma will be denoted by T max , which we will 
assume to be less than Asm- The reheating temperature, at which the thermal plasma domi- 
nates the energy density of the universe, will be denoted by Tr. In the simplest scenarios for 
reheating, we have the relation T^ ax ~ ^jT^Vj. Note that since tensor perturbations have not 
yet been observed in the cosmic microwave background, an upper bound has been placed on 
Vi of about (10 16 GeV) 4 . 

The vacuum decay rate induced by thermal fluctuations is again calculated by expanding 
around a bounce solution. The Higgs potential now involves a thermal contribution in addition 
to the one at T = 0: 

V tot (H) = ^p-H 4 + V TH (H; T). (32) 

The quartic potential with negative A does not have an SO(3)-symmetric bounce solution. 
A bounce exists only for the total potential fl32l) . The potential does not have a conformal 
symmetry at all, and as a result, bounces to field values of order T dominate in the decay rate, 
and the bounce action depends only on A(T). This is to be contrasted with the situation in 
section 2, in which tunneling to all scales had to be summed up, and the entire RG trajectory 
A(/z) was important for tunneling at any given time. It should be noted, though, that since the 
temperature of the universe changes, the full RG behavior of up to fi = T max will still be 
important in this section, through the time dependence of the temperature. 

We will now perform an analysis at a somewhat low level of approximation to get an overall 
qualitative picture of what is going on. For precise quantitative results, numerical calculations 
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will be used, as can be seen in what follows. We will employ the approximation adopted by 
[23] ; for g e &H <C T, the high-temperature (H/T) expansion of the potential may be used. The 
leading term of the 1-loop thermal potential is then the thermal mass term 

V® * \gl^H\ (33) 

where 

flSr = y 2 (j9 2 y + \ 9 l + 3/> 2 ) • (34) 

For A(T) negative and <C |A(T)|, the potential barrier for the transition occurs at around 
H ~ g c gT/ a/|A(T)|. Thus g e f[H is in this case much less than T, so that the high-temperature 
expansion is a good approximation. Reference [23] then obtained 

S3 6.0157rg eff (r) 

T - |A(T)| ' [6b) 
and a corresponding thermal decay rate per unit volume of 

/ q \ 3 / 2 s 

r[A(T);T]^T 4 ^J e~+. (36) 

In order to determine the probability of a point remaining in the false vacuum at a time t, 
one must consider the possibility of bubble nucleations throughout the point's past light-cone. 
The result for the fraction of the universe in the false vacuum at time t is then 



/(Ao; t) = exp (- ^ dh rlACTfa)); T(ty)]a{t x f^ If \ j , (37) 



where ti n it is the time right after the end of inflation, and a(t) is the scale factor of the expansion 
of the universe. In fact, the same equation was used in section 2 as a precursor to ffTUl) . Here, 
however, T changes as the temperature of the universe changes, and so we must now use the more 
general expression. Assuming a standard thermal history after the end of inflation, namely, a 
matter-dominated era of inflaton coherent oscillations, followed by radiation domination and 
matter domination (and further dark-energy domination), we obtain 

f(X ;t) ~ exp - / _r[A(T);T] 



t r T T 4 \T J \T E 

Tr dT rw , , Mpi {T(t)\ 3 r 

— r A T ; T — ^ \-zr-) t 3 \, (31 

A(T)<o T T 2 \ T J 
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assuming that t is in the recent matter-domination era0 Here, T eq is the temperature of matter- 
radiation equality, and T(t) is the photon temperature at time t. We have used T 4 oc l/a(t) 3 / 2 
and dt ~ — {MpiT R [/T i ){dT/T) during the inflaton dominated era [25J. Coefficients of order 
unity are ignored in the above expression for /, for the same reason as in section 2. 
One can introduce a critical value of the Higgs quartic coupling A C (T) by 



r[A c (T);T]r 3 ^ 
Mm 

J^2 



Tr 
T 



10 



T[X C (T)-T]r- 



Tr 

n 

T 



T R <T, 



T < Ti 



R, 



(39) 
(40) 



whereof T ~ 2.73 °K ~ 2.35 x 10 4 eV and r = 10 10 yrs. By taking a logarithm and using (1361) . 
we hava 15 ! 



6.0vr^ cfr (T) , / Mpi 3 / S 3 



7 In 



T 



or equivalently A C (T) 



6.07T5( efr (T) 



243 - In £ 



7 hi: 



(41) 
(42) 



The term [71n(T/T R )] is in square brackets because it should be included only when T > T R . 
The critical line of stability A C (T) is shown in Figure El along with some RG trajectories of the 
Higgs quartic coupling for m t = 171 GeV. If an RG trajectory passes below the critical line at 
an energy scale fi < T max , we will have /(Ao; £) <C 1, so that most patches of the universe with 
the corresponding Higgs boson mass would have decayed to the true vacuum by now. 

In Eq. ( 142]) . g e g(T), as calculated from ( 1341) . lies in the range 0.3-0.4 for a wide range of T. 



At the same time, |A C (T)| ranges from 0.03-0.04, so that |A C (T)| is not much bigger than g 4 ff . 
For this reason, the high-temperature expansion of the thermal potential is not particularly 
good, and consequently, 6.07rg e ff/|A| is not actually a very good estimate for the bounce action 



13 If t is in the era of cosmological constant dominance that follows after the recent matter-dominance era, 
then T(t) 3 t 3 is replaced by T^Hq 3 for t t', where T' eq and t' eq are the temperature and epoch of matter- 
cosmological constant equality, and Hq is the Hubble constant of the cosmological constant. Note that the 
t-dependence completely disappears from /(Ao;i). 

14 Factors of order unity are not important here, since these equations balance a very small T against a very 
large t 3 Mpi/T 2 . In particular, it is not important whether we are in the dark-energy dominated era or not, or 
which values to use for r and To. A factor of 2-3 changes A C (T) by half percent, which corresponds to less than 
0.1 GeV change in the Higgs boson mass. 

15 The number "243" contains (3/2) ln(243/27r) ~ 5.5 that comes from the 1-loop functional determinant. 
Note that AS of section 2, also calculated from a 1-loop functional determinant, contributed a similar order of 
magnitude, ~ —10 [4]. 
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Figure 6: The critical line of stability due to vacuum decay by thermal fluctuations, shown 
as solid lines terminating at T max . Two cases are shown, both with T max = 10 13 GeV. The 
upper solid line, case (I), corresponds to instantaneous reheating with Tr = T max . Varying 
Tr = T max away from 10 13 GeV will simply change the termination point of the line. The lower 
solid line, case (L), corresponds to late reheating with T R = 10 10 GeV. Varying T R will alter 
the position of the bend in the line, while altering T max will alter the termination point. The 
two lines coincide for T < 10 10 GeV. The long dashed line near the bottom is A C (M) discussed 
in section 2. Four RG trajectories of A(/i) are also shown, as dotted lines, with Higgs masses 
of (106.3, 113.1, 113.9, 116.5) GeV. The lower panel is a blow-up of the upper panel. The low 
temperature value of the critical Higgs mass, 106.3 GeV, is raised to 113.9 GeV in case (L) and 
116.5 GeV in case (I). 9 „ 



Sz/T [23]. The critical line of stability in Figure El was therefore obtained from a numerical 
calculation, the details of which are explained in the appendix. On the other hand, expression 
(|42|) does capture the qualitative behavior of A C (T) shown in the figure. For T < T R , the 
asymptotically free behavior of the SU(2)^ and top Yukawa couplings are more important than 
the —(1/243) ln(T/t>) term in the denominator, and |A C (T)| slowly decreases as T increases. For 
T > Tr, however, the extra —(7/243) ln(T/T^) term in the denominator is more important, 
and |A C (T)| increases. 

Figure shows A C (T) for two scenarios: one for T max = Tr = 10 13 GeV, and the other for 
Tr = 10 10 GeV and T max = 10 13 GeV. In the instantaneous reheating scenario at T = 10 13 
GeV (I), an RG trajectory of A(/i) with m# = 116.5 GeV (upper dotted curve in Figure [6b) 
touches the line of A C /} (T) at T ~ 10 13 GeV. Patches of the universe with Hide's boson masses 
smaller than 116.5 GeV would have mostly decayed because their RG trajectories pass below the 
A C 7) (T) line at some temperature less than T max . In the late reheating scenario with Tr = 10 
GeV (L), the RG trajectory for m H = 113.9 GeV (middle dotted curve in Figure Eb) touches 
Ac (T) at a temperature in between T R and T max . Thus, in this case patches of the universe 
with mn < H3.9 GeV typically decayed to the true vacuum within the inflaton dominated era. 
It is clear from Figure [6^ that the constraint from vacuum decay due to quantum fluctuations, 
shown as a long dashed line, is of negligible importance for both (I) and (L) scenarios discussed 
above, as well as many other cases with high Tr. 

4.2 The High Temperature Prediction for the Higgs Boson Mass 

When the a priori distribution -P(Ao) of the Higgs quartic coupling is highly peaked toward 
negative values, a balance between /(A ; t) and -P(A ) determines the most likely value of the 
Higgs boson mass. Since the double-exponential cut-off in / is so sharp, again the point of 
maximum probability occurs essentially at the cut-off, as in section 2. Thus, the phenomeno- 
logical lower bound on the Higgs boson mass may again become the prediction of a landscape 
scenario. 

Let us now consider a case with Tr = 10 10 GeV as an example. If T max is also 10 10 GeV, 
then / sets a cut-off A(10 10 GeV) > A c (10 10 GeV) ~ -0.026, as we see from Figure Eb. The 
RG trajectory with A(/i = 10 10 GeV) = —0.026 becomes the landscape prediction, which corre- 
sponds to ~m~H = 113.1 GeV. If T max is a little higher than 10 10 GeV, patches with tuh = 113.1 
GeV have mostly decayed during the inflaton-oscillation dominated era, because the RG tra- 
jectory lies below the critical line of stability \[ L \t). Thus the prediction for tuh goes up 
as T max increases. For sufficiently high T max , however, this sensitivity to T max disappears; the 
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RG trajectory with m# = 113.9 GeV touches \i L \T) at around 10 12 GeV. As long as T max is 
higher than this temperature, the landscape prediction will remain at ff%H = 113.9 GeV. Figure 
[7] shows the prediction for the most probable Higgs mass % in these two limits, the instanta- 
neous reheating scenario, with T max = Tr, and the high T max ^> Tr scenario. The prediction 
will always lie in between these two lines for any T max . 
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Figure 7: The prediction for the central value of the Higgs boson mass as a function of the 
reheating temperature Pr. The upper line is for the T max 3> Pr scenario, and the lower line for 
the instantaneous reheating scenario, T max = Tr. We assume that the top Yukawa coupling is 
fixed in the landscape, so that m t = 171 GeV. Note that the ±6 GeV error bar in (fT5l) . coming 
from threshold corrections at the weak scale etc., is not shown. 



Now, for sufficiently high reheating temperature, Pr ^ 10 13 GeV, the prediction does not 
depend on T max at all. This is because the RG trajectories for A(/z) run almost horizontally at 
high energies, and touch the critical line of stability A C (T) at the bend at T = Tr (see Figure 
ED- Vacuum decay becomes most effective at the epoch of reheating in this case, and it does 
not matter how high T max is. The fact that the RG trajectories for A run almost horizontally 
at high energies also explains why the prediction for m# levels off and depends only weakly 
even on Tr at high reheating temperatures. The Higgs quartic coupling does not run much at 
high energies because of the smaller top Yukawa coupling. 

Figure [6^ shows that the RG trajectory corresponding to the prediction of section 2 (the 
lowest trajectory) crosses A C (T) at around /i = 10 s GeV. Thus, the prediction of section 2 
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is unaffected if Tr = T max ^ 10 8 GeV in the instantaneous reheating scenario, as shown by 
point P2 in Figure UJ The highest reheating temperature for which the prediction of section 2 
persists actually depends on the relation between T R and T max . The higher T max is, the lower 
the maximum value of Tr for which the vacuum prediction applies. In the limit T max 3> Tr, 
this value of T R drops to around 10 5 GeV, as shown by point P\ in Figure [7J 

Finally, we have a few remarks about the uncertainties associated with the calculations of 
this section. The landscape prediction f[T5"j) in section 2, was derived from a value of the Higgs 
quartic coupling at a very high energy scale, M dom . The three sources of uncertainties quoted 
in ([TBI) , namely, the measured top quark mass, QCD coupling and calculation of electroweak 
threshold corrections, potentially changed the relation between \(Md om ) arid m#. The predic- 
tion in Figure [7J from the thermal scenario came from couplings renormalized at a lower energy 
scale ~ T Md om - However, the running of A(/x) takes place mainly at low energy, where 
the top Yukawa coupling is large. Thus, the prediction for the Higgs mass at high reheating 
temperature is also subject to the same three uncertainties quoted in (TH)]) with almost equal 
sizes. The influence of the details of the 1-loop functional determinant was at most 1 GeV [1] 
for the prediction in f[T5"j) . and it will remain at that order of magnitude in Figure [7J as well, 
because the estimates of the 1-loop functional determinant are much the same; see footnote [T5l 

It should be noted that the prediction in the case of high reheating temperature is associated 
with an additional uncertainty: the precision of the calculation of the thermal effective potential. 
This issue is briefly discussed in the appendix. 

4.3 The Width of the High Temperature Prediction 

The width of the Higgs mass prediction for this section may be obtained similarly to that of 
section 2. We consider the total probability V = Pf, and again we may define A_ and A + to 
be the values of the quartic coupling at Asm for which V has fallen from its peak value by a 
factor of 1/e. Since the cut-off from / is so sharp, A_ is again a negligible distance from the 
peak. As in section 2, A + is shifted above the critical value A c by an amount |A c |/p, where here 
we are using the value of the critical quartic coupling run up to the scale Asm- P is defined as 



in section 2 by p = 



. This yields a possible upwards shift in the Higgs mass of 

Ao=A c 

550|A c |GeV ,„ , 

Am H pa . (43) 

p 

Here, |A C | ranges from about .02 to .05 depending on the particular thermal history after 
inflation. This gives an upwards shift between (~ 10GeV)/p and (~ 25GeV)/p. 
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4.4 Discussion 

In the case that SM patches of the universe were reheated to very high temperatures after 
inflation, our prediction for the central value of the Higgs boson mass is shown in Figure El 
Despite a very wide range in the cosmologies considered, the central value remains in the region 
(106-118) GeV. Furthermore, for Tr < 10 5 GeV thermal fluctuations are irrelevant and the 106 
GeV result from quantum tunneling applies, while for Tr > 10 12 GeV, the prediction is in the 
narrow range (117 ± 1) GeV. The uncertainties from the experimental values of m t and a s and 
from higher loop orders are closely similar to those in the non-thermal case. The peaking of 
the a priori distribution for Ao yields an upper width on the prediction that is higher than the 
central value by an amount between ~ 10 GeV/p and ~ 25 GeV/p, depending on the reheat 
temperature. 

The top mass prediction of section 3 is very significant. If the top Yukawa coupling ho 
is allowed to scan, does this prediction survive thermal fluctuations at high temperatures? If 
Tr > 10 8 GeV, then a large p will causes h to rise significantly above the observed value, since 
a larger h gives both a steeper quartic trajectory and a lower A C (T) curve. Thus, to maintain 
the highly significant m t prediction from section 3, we are led to a preference for low Tr. 

5 The Scanning Standard Model 

Consider the SM, minimally augmented with right-handed neutrinos to allow for both neutrino 
masses and leptogenesis, with all parameters scanning — the scanning SM. It is possible that 
hospitable parts of the landscape exist where parameters take on values that are very differ- 
ent from the ones we observe. This is well-illustrated by the weakless universe |26j . where 
(H) /Mpi m 1 while the Yukawa couplings for the light quarks and the electron, h u ^, s ,e, are all 
less than 10 -20 . We assume that, if such regions exist, they are less probable than our own uni- 
verse. For example, the weakless universe may be disfavored because of a very small probability 
to have four Yukawa couplings that are so small. In the landscape, the origin for v -C Mpi 
may be that a small value for v is the most cost effective way of keeping u, d, s and e light. 
Hence the question becomes: in our neighborhood of the landscape, how much of the parameter 
space of the scanning SM can be understood from environmental selection? We have argued 
that the entire Higgs potential is strongly selected, and now we turn to the Yukawa coupling 
matrices, where most of the parameters lie. In sections 2, 3, and 4 we have been exclusively 
focused on precise predictions for the Higgs and top masses. In this section we explore possible 
consequences of embedding our previous results into the scanning SM. Further assumptions are 
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Figure 8: The Yukawa couplings for the quarks and charged leptons scaled up to 10 18 GeV, 
and plotted logarithmically. The environmentally important e, u, d and t couplings seem to be 
outliers, while the remaining five couplings are all clustered in a two order of magnitude range. 



found to lead to intriguing order of magnitude estimates for a variety of quantities. 
5.1 Quark and Charged Lepton Masses 

Let us first ignore flavor violation and consider scanning of the nine charged fermion eigenval- 
ues. In Figure [8] we plot the nine Yukawa couplings of the SM at the renormalization scale 
Asm = 10 18 GeV. They fall into three well separated groups: the e, u and d Yukawas, selected 
to be small by the requirement that atoms exist [HT], the t coupling, selected to be large to 
avoid vacuum instability, and the remaining five Yukawa couplings that are environmentally 
irrelevant, and therefore reflect the underlying probability distribution on the landscape. It is 
striking that while the nine couplings span 5 orders of magnitude, the five irrelevant couplings 
span only 2 orders of magnitude. It therefore seems possible to make progress on the fermion 
mass hierarchy problem by assuming a universal fermion Yukawa eigenvalue probability distri- 
bution that is peaked around 10~ 3 , with a width of one to two orders of magnitude. Possible 
simple analytic forms for the distribution include 

P(h) =Ae v^s^j (44) 

with h = /io/lCT 3 and a ~ 0.75, and 

P(h) = A~he-~ h . (45) 

Thus the five irrelevant Yukawa couplings are scattered logarithmically about 10 -3 , while the 
four outliers are environmentally selected to be far from the peak of the distribution^ This 

16 With these distributions, the weakless universe is indeed much less probable than the observed SM. 



28 



interpretation of the data of Figure [S] seems so plausible, that we investigate further the role of 
electroweak phase stability in forcing the top coupling to large values. 

From Figure [2] it is clear that if a strongly peaked distribution -P(A ) pushes A towards its 
critical value A U|C , then vacuum instability is only avoided if there is a Yukawa coupling that is 
sufficiently large to RG scale A to positive values by the weak scale. This can be accomplished 
by a large top coupling via the term — Qhf of (]7j). However, with all Yukawa couplings scanning, 
the — term could be replaced by other large Yukawa contributions. One easily sees that it 
is more probable to have a single large Yukawa rather than two or more large Yukawas making 
significant contributions to the quartic RG equations. Furthermore, it is not sufficient for the 
single large Yukawa to be that of a lepton: it is only for a quark that QCD increases the 
Yukawa in the infrared, thus magnifying the effect on the scalar quartic coupling. The critical 
value for a lepton Yukawa coupling is slightly larger than for a quark, but this is insufficient 
to allow a critical quartic coupling at Asm to become positive by the weak scale. Hence, our 
simple interpretation of Figure [HI in terms of a universal P(h ) and environmental selection 
of outliers, predicts not only that there will be a single heavy fermion, but that it must be a 
quark. However, the choice of up-like versus down-like is random. Crucially, as we have shown 
in section 3, a precise and successful prediction for the top quark results] 17 ! 

If this picture of fermion masses is correct, then there is a further consequence for the u, d 
and e masses. Since they are selected to be quite far from their most probable values in the 
landscape, our patch of the universe is quite rare. This means that they are expected to be 
quite close to the maximum values consistent with environmental selection for the existence of 
complex atoms. Of course, this is sensitive to the actual probability distribution for the Yukawa 
couplings. For example, using the distribution (J32J) patches in the multiverse with the electron 
Yukawa double its observed value are more than an order of magnitude more probable than in 
our patch. On the other hand, for the distribution (H51) they are only twice as probable. 

It is not clear how the small CKM mixing angles could arise from environmental selection. 
Vacuum stability has selected an up-type quark to be heavy, and nuclear physics has selected 
a down- type quark to be light. In nature these two quarks are almost precisely in different 
weak doublets, so that V t d — 10~ 3 , but there does not appear to be a selection effect to explain 
this. This is similar, but numerically more severe, to the puzzle of why environmental selection 
would make V u d close to unity. Quark flavor violation appears to be governed, at least to 

17 Note that the distribution ([44]) gives a peaking parameter q ~ 2. While this has the right sign to limit 
Amt+, the size of A.mt+ from (|30|) seems too large. Hence, either P(h) drops more rapidly at large h compared 
to (T44|) . or Amt+ is made small by a large tunneling rate in the full theory above Asm- 
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some degree, by symmetries. Numerical simulations of Yukawa matrices have been made for 
very flat landscape probability distributions [27] . It would be interesting to perform numerical 
simulations using more peaked distributions, such as (1441) or (1451) . and incorporating elements 
of approximate flavor symmetries. 

5.2 Neutrino Masses, Leptogenesis and Inflation 

Since no large neutrino mass ratios or very small neutrino mixing angles have been measured, 
it is very plausible that the entries of the neutrino mass matrices, both Dirac and Majorana, 
are distributed according to some universal probability distribution of the landscape. Indeed, 
Neutrino Anarchy demonstrated in some detail how randomly generated neutrino mass matrices 
could account well for the observed masses and mixings [28]. The view of the charged fermion 
masses presented above has an important consequence for neutrino physics. We assume that the 
Dirac neutrino Yukawa couplings are governed by the same landscape probability distribution as 
the quarks and charged leptons. If there is little environmental selection acting on the neutrinos, 
then the typical Dirac Yukawa eigenvalue h v is expected to be of order 10~ 3 in magnitudecj 
Since no large hierarchy of neutrino masses is observed, it is reasonable to also assume that the 
three Majorana masses of the right-handed neutrinos are governed by a universal probability 
distribution, so that they are all expected to have a common order of magnitude, Mr. With 
h u ~ 10 -3 , there is a relation between Mr and the order of magnitude of the light neutrinos, 
fn u : 

m v 3 x 10 7 GeV / 7L 



eV M R \10~* 

In our patch of the multiverse, m v has been observed to be of order 10~ 2 eV, so that we predict 

M R « 3 x 10 9 GeV (47) 

in our patch. 

Leptogenesis is the only possible source of a baryon asymmetry in the scanning SM. Is it 
consistent with our framework for fermion masses? The lepton asymmetry is given by [29] 

Yl « >A (48) 
9* 

where ~ 100 counts the number of states in the thermal bath. The CP asymmetry e\ and 
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Sufficiently small that these couplings do not upset the prediction for the Higgs boson and top quark masses. 
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the washout factor k are roughly of order 

1d7t nf, Mpi 

respectively, where we have assumed comparable masses for the right-handed neutrinos. We 
are now assuming that all relevant entries of the neutrino Yukawa coupling matrix are of order 
h u . We see that the observed baryon asymmetry is roughly reproduced after using (147|) : 

n«io-"f g » Wo-"f^)7— V (so) 

V3xl0 9 GeVy V 10 / \ / 

This can be regarded as an indication that the Dirac Yukawa couplings of neutrinos are also 
subject to the universal distribution P(ho) of Yukawa couplings. 

The baryon asymmetry is roughly derived from the observed low-energy neutrino masses, or 
vice versa; either one of them can be used as an input so that the other is determined. Now the 
next question is whether we can understand them all together environmentally. Suppose that 
the Majorana masses of the right-handed neutrinos are scanned. The distribution of Mr may be 
a featureless flat distribution, like that of the cosmological constant; or, the baryon asymmetry, 
being proportional to Mr (as in (j50p ). may put more weight on larger values of Mr. We assume 
that the distribution of Mr is mildly peaked toward larger values. If the reheating temperature 
Tr is also scanned mildly, but with an a priori distribution peaking weaker than that of Mr, 
then both Mr and Tr would go hand in hand, pushed upward, because the baryon asymmetry 
is at least power-suppressed when Tr is much lower than Mr. Both Mr and Tr go up until Tr 
and T max are so high that the Higgs vacuum becomes unstable. If T max is also scanned, with 
an even milder a priori distribution, Tr will be pushed upward along the line from Pi to P 2 
in Figure [7] until it reaches the special point P 2 . For a large p, on the other hand, a reheating 
temperature higher than that is less probable, because a large p strongly favors a smaller Ao, 
which is environmentally allowed only in patches with lower Tr. Thus environmental selection 
for the desired electroweak vacuum, on the combined a priori probability distribution, leads to 
the special point P 2 as the most likely one to be observed, and hence to the predictions 

Tr w T max w 10 8 GeV. (51) 

The reheat temperature is then low enough that the predictions for the Higgs and Top masses 
given in sections 3 and 4 remain valid. 

Since T max pa ^/T^Vi, we predict the vacuum energy during inflation 

Vi Pa (10 8 GeV) 4 (52) 
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corresponding to a Hubble parameter during inflation of order 10~ 3 GeV. These speculations 
will clearly be disproved if a tensor perturbation is seen in the cosmic microwave radiation. 
From the reheat temperature we predict the decay rate of the inflaton 

r 7 ^10" 3 GeV. (53) 

Now that the rough scale of Mr is obtained without an observational input, low-energy neutrino 
masses of order 10 -2 eV and the baryon asymmetry are also purely predictions of environmental 
selection. 

Since in thermal leptogenesis the lightest right-handed neutrino mass can be slightly higher 
than the reheating temperature Tr, there is no conflict between Tr in fl5Tl) and Mr in (jUJ). Also, 
the neutrino Yukawa couplings contributing to the CP asymmetry e± may fluctuate bigger, and 
those appearing in the denominator of the washout factor k smaller, making the entire baryon 
asymmetry bigger. Thus an apparent discrepancy of one order of magnitude is nothing to 
worry about. Although the reheat temperature (1511) does not satisfy the phenomenological 
limit 2 x 10 9 GeV < Tr derived from the observed baryon asymmetry [29 J, our estimates 
are only order of magnitude ones. Furthermore, this limit does not apply when the right- 
handed neutrinos have highly degenerate Majorana masses. An environmental factor oc Yr 
may motivate such highly degenerate right-handed neutrinos. 

5.3 Gauge Couplings 

If all three SM gauge couplings scan, then a and a s are selected, while the weak mixing angle 
sin 2 9 is apparently environmentally irrelevant. However, the condition for the critical value for 
the top quark coupling, equation fl27|) . depends on the weak mixing angle. Hence it is possible 
that avoiding vacuum instability could select for the weak mixing angle in addition to the top 
quark and Higgs masses! If the a priori probability distribution has a stronger dependence on 
h t than sin 2 9 then we find the incorrect prediction of sin 2 9 ~ 0.38 at the weak scale. Hence it 
must be that the weak mixing angle is essentially determined by the landscape distribution. 

Gauge coupling unification is possible in the scanning SM, with Asm near Mpi, if threshold 
corrections are of order 10-20%. Since the threshold corrections to the gauge coupling constants 
arise from dimension-5 operators, while the leading order correction to the Higgs potential ([6]) 
is a dimension-6 operator, corrections to gauge coupling unification can be as large as 10-20% 
without losing the 1%-level precision of the predictions in earlier sections. 
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5.4 Dark Matter 



The final parameter of the SM that could scan is the strong CP parameter of QCD. It is 
attractive to promote this to the axion field, with an associated symmetry breaking scale above 
Asm- It has been argued that scanning the primordial value of this axion field not only allows 
axionic dark matter, but could explain the approximate equality of dark matter and baryon 
energy densities [33] . 

The scanning SM is frustrating in the sense that there appears little that is amenable to 
experimentation beyond the Higgs mass prediction. However, the question of flavor in both 
quark and neutrino sectors is left unresolved, and in the neutrino sector more thought needs to 
be given to the environmental selection effects arising from leptogenesis. 

We caution the reader that there are potential difficulties with the scanning SM that we do 
not address in this paper. For example, particle physics of the SM is completely unchanged if 
m 2 and the QCD scale Aq CD are increased by the same factor, and corresponding changes are 
made to the dimensionless couplings so that their values at the scale m? are kept fixed. The 
environmental selection of the weak scale to be as small as 100 GeV would then be undermined, 
unless some probability distribution can be found to counteract the one favoring large m 2 , or 
some other environmental selection involving the Planck scale, such as big bang nucleosynethe- 
sis, limits such a scanning direction. Similar questions arise for the environmental selection of 
the cosmological constant when there is simultaneous scanning of other cosmological quantities, 
such as the density perturbations. This has been addressed by a consideration of other selection 
effects [30]. 

6 Conclusions and Discussion 

There are many mass scales that play an important role in physics, for example the electron 
and proton masses, but from a fundamental viewpoint two are key: the scale of electroweak 
symmetry breaking, (H) = v, that sets the mass scale of the quarks and leptons as well as the 
W and Z bosons, and the Planck mass, Mpi, that describes the gravitational coupling and is the 
largest scale that enters physics. We take the viewpoint that the Planck scale is fundamental 
and the weak scale is somehow to be derived, although the alternative is also possible, so that 
the key question becomes the origin of the small dimensionless parameter v/Mpi. 

Over the last 30 years, many frameworks have been developed that can explain this small 
ratio by introducing various symmetries. The LHC is expected to play a crucial role in dis- 
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tinguishing between these frameworks, for example between supersymmetry and new strong 
forces at the TeV scale. The alternative viewpoint, that v/Mpi is environmentally selected 
[8] and has nothing to do with extended symmetries, has received relatively little attention, 
perhaps because it leads to no new physics, beyond a light Higgs boson, to be discovered at 
LHC. Of course it is possible that in this case the LHC might discover physics unrelated to 
electroweak symmetry breaking, perhaps associated with either dark matter or coupling con- 
stant unification, but it is hardly to be expected. Dark matter is not necessarily related to the 
electroweak scale; axion dark matter with an environmentally chosen relic density is an exam- 
ple. Gauge coupling unification requires new particles other than a Higgs doublet, but those 
particles may be around the Planck scale, for instance. If the origin of v/Mpi is environmental, 
then arguments for physics discoveries beyond the SM at LHC become tenuous. 

Another fundamental mass scale of nature is the cosmological constant, A 4 . To date, theories 
beyond the SM have failed to provide a symmetry understanding for why A/Mpi is so small: 
why it is so close to zero, and what determines the order of magnitude of the deviation from zero. 
A symmetry might assign a special meaning to a vanishing cosmological constant, and it might 
be possible to understand the observed size of the dark energy by constructing theories that 
implement a seesaw relation A m v 2 /Mpi. On the other hand, environmental selection explains 
both the extreme smallness of A and its order of magnitude, as measured by the dark energy 
density. The dark energy and the ratio v/Mpi are environmentally selected if two parameters 
of the renormalizable Higgs quartic potential are scanned. These successes motivate us to think 
of scanning the entire Higgs potential. In this case, the Higgs quartic coupling, A, and therefore 
the physical Higgs mass, scans. By deriving a prediction for the Higgs boson mass, the idea of 
the scanning Higgs potential can be tested. 

Environmental selection requires an "edge" [2T], [32] — a surface in the space of scanning 
parameters, such that on one side of the surface the formation of desired complex structures is 
greatly suppressed. We do not know enough about conceivable life forms to derive the precise 
types of complexity to be selected. However, the more general the requirement, the more 
plausible selection becomes. The selection of small values for A/Mpi results from an edge that 
corresponds to formation of non-linear structures late in the universe. As A/Mpi is increased 
beyond this edge, the probability of such structures forming is suppressed because they are 
subjected to inflation before they go non-linear. On the other hand, a small value for v/Mpi is 
selected by the requirement that atoms exist. The requirement that atoms exist is much more 
specific than the requirement that large non-linear structures form, but it is not unreasonable. 
For our prediction of the Higgs boson mass, the relevant edge is the phase boundary that 
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separates a metastable phase with small (H) /Mpi = v/Mpi from a phase with (H) /Mp\ ~ 1. 
This is less specific than requiring the existence of atoms: life probably requires objects that 
contain a large amount of information, and this may be accomplished best in a phase with a 
large Mpi/v ratio. 

An edge that allows environmental selection is not sufficient to make predictions; some 
knowledge of the a priori probability distribution for the scanning parameters is needed. In 
the absence of a calculation from a fundamental theory of the landscape, this requires an 
assumption. For the cosmological constant and weak scales, this assumption is extremely 
mild; it is sufficient to assume that the distributions are roughly flat and featureless, so that 
small values of these scales are exceedingly rare. For our Higgs mass prediction a non-trivial 
assumption is necessary: the probability distribution P(\) must be sufficiently peaked at low 
coupling that A is expected to be near the metastability boundary. 

Since our edge is a phase boundary, the position of the boundary, and therefore the Higgs 
mass prediction, depends on the thermal history of the universe, in particular on the reheat and 
maximum temperatures after inflation, Tr and T max , two of the most important parameters of 
post inflationary cosmology that are still largely unconstrained, and consequently one might 
expect that little can be said about the Higgs mass. In fact, no matter what the values of these 
two temperatures the most probable value of the Higgs mass is raised from 106 GeV to at most 
118 GeV. Furthermore, if T max = Tr & 10 8 GeV the prediction is very insensitive to cosmology, 
as the danger of a phase transition is highest today and was negligible in the early universe, 
giving a central value 



m, = 10. GeV + 6 G eV (^^) - 2.6 GeV («^) ± 6 GeV. 



(54) 



with an upward statistical fluctuation of (25/p) Gev. For Tr > 10 12 GeV, nucleation at high 
temperature dominates, giving m# = 117 ± 1 GeV for any value of T max , with an upper fluctua- 
tion between ~ (10/p) GeV and ~ (25/p) GeV from the a priori distribution, depending on the 
reheat temperature. The full dependence of the prediction on Tr and T max is shown in Figured 
The parameter p is the logarithmic derivative of P(X) evaluated at the phase boundary. Further 
uncertainties, displayed in equations (1151) and (fT8"|) . arise from the experimental uncertainty in 
the QCD coupling and the top quark mass, and also from use of perturbation theory in RG 
scaling of couplings and in relating couplings to physical masses. Ultimately, these uncertain- 
ties can be reduced; but the width of the Higgs boson mass distribution from the multiverse, 
parameterized by p, cannot be reduced. This width is comparable to the present uncertainties 
if p pa 3, while if p > 10 it is practically negligible compared with other uncertainties. 
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The top quark plays a crucial role in electroweak symmetry breaking via the RG evolution 
of the quartic coupling A. If the top Yukawa coupling also scans, a further mild assumption 
also allows a prediction for the top quark mass. For Tr < 10 8 GeV, the most probable top 
quark mass is 

m t = [176.2 ±3 + 2.2 log 10 (A 5A/ /10 18 GeV)] GeV. (55) 

To limit fluctuations above this requires a further assumption, while downward fluctuations are 
35 GeV/ v /p. The striking success of this result suggests that environmental selection may be 
at work in the Higgs potential, and that we may be lucky enough for p to be sufficiently large 
to give a precise Higgs mass prediction. Maintaining this prediction in the presence of thermal 
fluctuations also implies that Tr ^ 10 8 GeV. 

We have argued that these Higgs and top mass predictions could occur in the scanning SM, 
where all SM parameters scan. Indeed a zeroth-order understanding of the charged fermion 
masses follows if all Yukawa couplings have a universal probability distribution centered on 
about 10~ 3 . An interesting extension to neutrino masses and leptogenesis follows if right- 
handed neutrino masses scan, with a preference for larger values, and if Tr and T max scan 
with mild distributions. The broad order of magnitude of the light neutrino masses and the 
baryon asymmetry are correctly predicted, while the right-handed neutrino masses, the reheat 
temperature and the maximum temperature are all predicted to be of order 10 8 -10 9 GeV. 

The possibility of precise predictions for the Higgs and top masses by environmental selec- 
tion from a landscape illustrates that it may be possible to do physics without symmetries. 
Theories based on symmetries yield predictions because the symmetries limit the number of 
free parameters. One chooses a highly symmetric model and studies the resulting predictions. 
On the landscape two ingredients are needed for a "bottom-up" prediction: an edge that allows 
environmental selection, and an assumption that the a priori probability distribution is pushing 
parameters towards this edge. Predictions follow, quite literally, from living on the edge. Thus 
inventing models is replaced by finding edges and inventing probability distributions. There 
are two difficulties with this approach: relevant edges are hard to come by, and even if environ- 
mental selection is occurring at an edge, we can only discover it if the probability distribution 
happens to be favorable. 

If the LHC discovers a light Higgs boson in our predicted range, and no sign of any physics 
beyond the SM, then a precise numerical test of environmental selection could follow from 
further developments in both theory and experiment. This will require a reduction in the 
experimental error bars of the QCD coupling and the top quark mass, and improved accuracy 
of calculations of both RG scaling and extracting pole masses from running couplings. In the 
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case that the phase boundary is determined by thermal fluctuations at high temperatures, a 
refined calculation of the temperature dependent effective potential will be needed, in particular 
including the effects of Higgs quanta in the thermal bath. 



We conclude by addressing those who are sceptical of physics arguments from environmental 
selection. We have shown how precise predictions for the Higgs and top masses arise from 
assumptions about the landscape. If these precise predictions are found to be false, then our 
theoretical construct will be experimentally disproved. If they are found to be correct, then the 
SM is indeed on the edge of electroweak vacuum instability. It would then be up to you to find 
an alternative, more convincing theoretical framework to explain this fact. In this situation, we 
suspect that effort will be focussed on extracting further predictions from the landscape, and 
on understanding how the landscape might arise from string theory. 
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A A Little More about The Thermal Effective Potential 

At the 1-loop level, the effective potential in a thermal bath V to t{H) can be split into a sum of 
the effective potential at T = and a thermal 1-loop contribution, Vth{H;T). The thermal 
contribution is given by [221 1231 123] 



V TH (H;T) 




where 




(57) 



(58) 
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Figure 9: Effective potential in the thermal bath V tot {H) for different values of A, (from left) — 1, 
—0.1 and —0.01. The potential with (156]) is drawn by solid lines, and with its approximation 
(155]) by dashed lines. Using g L = 0.570, g Y = 0.402 and h = 0.547, we find g cS ~ 0.382 and 
g*s ~ 0.02 for T = 10 10 GeV and m H = 120 GeV. 



Here, only contributions of loops of W-bosons, Z-bosons and top quarks are taken into account, 
and not of the Higgs field itself. For Higgs field values much less than the temperature, this 
thermal contribution to the effective potential is approximated by a free energy piece plus 
quadratic term. On the other hand, for T <^ H, Vth{H\ T) in fl56l) is exponentially small. 

For a negative Higgs quartic potential with large |A(T)|, the potential barrier for the vacuum 
transition is at a small field value of H . Hence, the overall shape of the potential is almost 
properly obtained even if the thermal contribution is approximated by the quadratic term 
(see Figure [Uti). On the other hand, for a small |A(T)|, the potential barrier is outside the 
range where the high-temperature expansion is valid, and Vtot calculated using the quadratic 
approximation for the thermal contribution (133]) is quite different from the right shape of the 
potential (see Figure [9b). The quadratic approximation V^(H;T) in (133]) continues to grow 
for large H, but the true form of Vth(H; T) levels off for H T. Thus, the potential using the 
quadratic approximation tends to overestimate the height and width of the potential barrier 
of the vacuum transition for small |A|. From Figure O a-c, the quadratic approximation of the 
thermal potential is good for |A] ~ 1 3> g^s ~ 10~ 2 , but clearly not for |A| ~ 10 -2 . 

Since the quadratic approximation of the thermal potential Vth{H;T) overestimates the 
height and width of the potential barrier for small |A|, the bounce action S3/T ~ 6.015 7rg e fj/|A| 
based on the approximation is also an overestimation. Figure [TUk shows the bounce action 
calculated numerically with the full form of the potential Vth(H; T) in ([56]) and with V^(H; T) 
in ([33]) . Figure [TOb shows the ratio between them. As expected from the difference in the shape 
of the potential in Figure El the bounce action is overestimated in (1351) — by 20-30% from Figure 
[TUb for |A| ~ 0.02-0.04, the range of practical interest. Thus, if the landscape prediction of the 
Higgs boson mass were calculated upon the quadratic approximation, A C (T) would have been 
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Figure 10: The bounce action S$/T with (dashed) and without (solid) quadratic approximation 
of the thermal potential. 



20-30% smaller, corresponding to A at the weak scale being 10% smaller, and the prediction 
for m H being 5% lower. It is not appropriate to use the quadratic approximation to calculate 
A C (T) and to obtain a precise landscape prediction. All the calculations in section 4 use the 
potential (151)1) . and the bounce solution and action were obtained and calculated numerically. 
We used m# = 120 GeV as the boundary condition of the 2-loop RG equations in order to 
calculate gy, Ql and h at high energy scale used in (1561) . but there is no impact on the result 
shown in Figured even if m# = 130 GeV is used instead. 

On the other hand, Figure [TOa also tells us that (1351) captures the qualitative aspects very 
well; the numerically calculated bounce action (solid line) behaves almost the same way as (1351) . 
Thus, the expression ( 1351) as an approximation of the bounce action is still quite useful when 
thinking of qualitative issues. 

The biggest problem with the treatment of the thermal potential so far is that Higgs loops 
have not been included in (156|) F^I The contribution to g^ s has also been omitted; it would have 
been 

9*. = ^ (j9y + \qI + 3^ 2 + 6a) . (59) 

It is technically quite involved to incorporate the Higgs loop contribution to the thermal poten- 
tial, when the high-temperature expansion is no longer valid. This is beyond the scope of this 
paper. Instead, let us try to get a feeling for the size of the possible effects of the Higgs loop 
contribution, using expressions in the quadratic approximation, which seems to be qualitatively 

19 The contribution from bottom quark loops is irrelevant to the calculation of the bounce action so long as 
^«|A|. 
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valid. The 6A contribution in g^ s leads to a 5% change in g e g_ for |A| = 0.03. A shift in A C (T) of 
roughly 5% maintains the value of the bounce action (l3~5l) . corresponding to a 2% change in the 
Higgs quartic coupling at the electroweak scale, and 1% change in the Higgs mass prediction. 
This estimate, however, heavily relies on the high-temperature expansion, and thus may not 
be particularly accurate. 
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